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ABSTRACT 

We have undertaken a program to observe emission lines of [S IV] 10.51, 
[Ne II] 12.81, [Ne III] 15.56, and [S III] 18.71 /im in a number of extragalactic 
H II regions with the Spitzer Space Telescope. Here we report our results for 
the nearly face-on spiral galaxy M83. A subsequent paper will present our data 
and analysis for another substantially face-on spiral galaxy M33. The nebulae 
selected cover a wide range of galactocentric radii (Rg). The observations were 
made with the Infrared Spectrograph in the short wavelength, high dispersion 
configuration. The above set of four lines is observed cospatially, thus permitting 
a reliable comparison of the fluxes. From the measured fluxes, we determine the 
ionic abundance ratios including Ne ++ /Ne + , S 3+ /S ++ , and S ++ /Ne + and find 
that there is a correlation of increasingly higher ionization with larger R^. By 
sampling the dominant ionization states of Ne and S for H II regions, we can 
approximate the Ne/S ratio by (Ne + + Ne ++ )/(S ++ + S 3+ ). Our findings of 
ratios that significantly exceed the benchmark Orion Nebula value, as well as 
a decrease in this ratio with increasing Kg, are more likely due to other effects 
than a true gradient in Ne/S. Two effects that will tend to lower these high 
estimates and to flatten the gradient are first, the method does not account for 
the presence of S + and second, S but not Ne is incorporated into grains. Both 
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Ne and S are primary elements produced in a-chain reactions, following C and O 
burning in stars, making their yields depend very little on the stellar metallicity. 
Thus, it is expected that Ne/S remains relatively constant throughout a galaxy. 
We stress that this type of observation and method of analysis does have the 
potential for accurate measurements of Ne/S, particularly for H II regions that 
have lower metallicity and higher ionization than those here, such as those in 
M33. Our observations may also be used to test the predicted ionizing spectral 
energy distribution (SED) of various stellar atmosphere models. We compare 
the ratio of fractional ionizations <Ne ++ >/<S ++ > and <Ne ++ >/<S 3+ > vs. 
<S 3+ >/<S ++ > with predictions made from our photoionization models using 
several of the state-of-the-art stellar atmosphere model grids. The overall best fit 
appears to be the nebular models using the supergiant stellar atmosphere models 
of Pauldrach et al. (2001) and Sternberg, Hoffmann, & Pauldrach (2003). This 
result is not sensitive to the electron density and temperature range expected 
for these M83 nebulae. Considerable computational effort has gone into the 
comparison between data and models, although not all parameter studies have 
yet been performed on an ultimate level (e.g., in the present paper the stellar 
atmosphere model abundances have been fixed to solar values). A future paper, 
with the benefit of more observational data, will continue these studies to further 
discriminate how the ionic ratios depend on the SED and the other nebular 
parameters. 

Subject headings: ISM: abundances, H II regions, stars: atmospheres, galaxies: 
individual (M83) 

1. Introduction 

Most observational studies of the chemical evolution of the universe rest on emission line 
objects, which define the mix of elemental abundances at advanced stages of evolution as well 
as the current state of the interstellar medium (ISM). Gaseous nebulae are laboratories for 
understanding physical processes in all emission-line sources and probes for stellar, galactic, 
and primordial nucleosynthesis. H II regions are also among the best tracers of recent star 
formation. 

The presence of radial (metal/H) abundance gradients in the plane of the Milky Way 
is well established in both gaseous nebulae and stars (e.g., Henry & Worthey 1999; Rolle- 
ston et al. 2000). Radial abundance gradients seem to be ubiquitous in spiral galaxies, 
though the degree varies depending on a given spiral's morphology and luminosity class. 
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The gradients are generally attributed to the radial dependence of star formation history 
and ISM mixing processes (e.g., Shields 2002). Thus, the observed gradients are a major 
tool for understanding galactic evolution (e.g., Hou et al. 2000; Chiappini et al. 2001; Chiap- 
pini et al. 2003). The premise is that star formation and chemical enrichment begins in the 
nuclear bulges of the galaxies and subsequently progresses outward into the disk, which has 
remained gas-rich. The higher molecular gas density in the inner regions produces a higher 
star formation rate, which results in a relatively greater return to the ISM of both "primary" 
a-elements (including O, Ne, and S) from massive star supernovae, and "secondary" elements 
like N. Secondary nitrogen is produced by CNO burning of already existing carbon and oxy- 
gen in intermediate-mass stars and is subsequently returned to the ISM through mass loss. 
However, because chemical evolution models have uncertain input parameters, and because 
details of the abundance variations of each element are uncertain, current understanding of 
the formation and evolution of galaxies suffers (e.g., Pagel 2001). 

Studies of H II regions in the Milky Way are hampered by interstellar extinction. For the 
most part, optical studies (e.g., Shaver et al. 1983) have been limited to those H II regions at 
galactocentric radius > 6 kpc (predicated on R & = 8 kpc) because H II regions are very 
concentrated to the Galactic plane. Here extinction becomes severe with increasing distance 
from Earth. Observations using far-infrared (FIR) emission lines have penetrated the Kg < 
6 kpc barrier. Surveys with the Kuiper Airborne Observatory (KAO) by Simpson et al. 
(1995), Afflerbach et al. (1997), and Rudolph et al. (2006) have observed 16 inner Galaxy 
H II regions. With the Infrared Space Observatory (ISO), Martin- Hernandez et al. (2002a) 
observed 13 inner Galaxy H II regions covering FIR and also mid-IR lines. A major finding 
of these studies is that inner Galaxy H II regions generally have lower excitation (ionization) 
compared to those at larger Kg- This holds for both heavy element ionic ratios ++ /S ++ 
(Simpson et al. 1995) and Ne ++ /Ne + (Simpson & Rubin 1990; Giveon et al. 2002), and also 
He + /H + measured from radio recombination lines (Churchwell et al. 1978; Thum, Mezger, 
& Pankonin 1980). Whether the observed increase in excitation with increasing R^ comes 
entirely from heavy element opacity effects in the H II regions and stellar atmospheres, or 
also from a gradient in the maximum stellar effective temperature, T e s, of the exciting stars 
is still a point of controversy (e.g., Giveon et al. 2002; Martin-Hernandez et al. 2002b; Smith, 
Norris, & Crowther 2002; Morisset et al. 2004). 

It has become clear that nebular plasma simulations with photoionization modeling 
codes are enormously sensitive to the ionizing spectral energy distribution (SED) that is 
input (e.g., Simpson et al. 2004, and references therein). These SEDs need to come from 
stellar atmosphere models. Stellar atmosphere modelers are paying increasing attention 
to the usefulness of nebular observations, particularly of H II regions, in validating and 
constraining their models (e.g., Stasihska & Schaerer 1997). Though it may be somewhat 
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chauvinistic to say so, much of this increased attention stems from work we did pointing 
out the "[Ne III] problem" and possible solutions (Simpson et al. 1995; Rubin et al. 1995; 
Sellmaier et al. 1996). To produce Ne ++ requires ionizing photons > 41 eV. At energies 
exceeding this ionization potential (IP), modern stellar atmosphere model SEDs are espe- 
cially divergent. A critical test of the validity of stellar atmosphere models of hot stars is 
whether H II region models produced with these atmospheres predict line fluxes that agree 
with observations. The "[Ne III] problem" is that the observed Ne ++ /0 ++ ratio significantly 
exceeds model predictions and remains relatively constant over a large range of H II region 
excitation as gauged by the ++ /S ++ ratio. The original observational basis is KAO FIR 
measurements of [Ne III] 36, [O III] 52, and [S III] 33 /j,m lines in Galactic H II regions (see 
figure 3 in Simpson et al. 1995). 

While Sellmaier et al. (1996) believed that they had solved the [Ne III] problem when 
they obtained a good fit to the FIR data by using non-LTE atmospheres with winds computed 
with Pauldrach's code as it then existed, we recently investigated the problem and concluded 
that it still exists (Simpson et al. 2004). Our H II region models using non-LTE stellar 
atmospheres with winds from Pauldrach, Hoffmann, & Lennon (2001) and from atmosphere 
models by both Smith et al. (2002) and Sternberg, Hoffmann, & Pauldrach (2003) using 
Pauldrach's WM-BASIC code predict much lower Ne ++ /0 ++ than observed for "Dwarf 
atmospheres with T eff < 40, 000 K and for "Supergiant" atmospheres with T eff < 35, 000 K. 
Moreover, the models without winds do a poorer job of reproducing the observations than 
those that include winds. 

In this paper, we develop new observational tests of and constraints on the ionizing 
SEDs that are predicted from various stellar atmosphere models. We do this by utilizing 
Spitzer Space Telescope (SST) observations of H II regions in the spiral galaxy M83. M83 
(NGC 5236) is one of the closest (distance 3.7 Mpc) and brightest spirals (SBc II), as well 
as being nearly face-on (i = 24°) (de Vaucouleurs et al. 1983). M83 particularly interested 
us because of its high metallicity (at least twice solar, e.g., Dufour et al. 1980; Bresolin & 
Kennicutt 2002). High metallicity should correlate with lower-ionization H II regions. From 
our experience with the [Ne III] problem, this is likely to provide the most stringent test 
matching theory with observations. With the Infrared Spectrograph (IRS) on the SST, we 
can observe four emission lines that probe the dominant ionization states of neon and sulfur 
in these H II regions. These lines are: [Ne II] 12.81, [Ne III] 15.56, [S III] 18.71, and [S IV] 
10.51 /im. Crucially, SST permits the simultaneous, cospatial observation of these four 
lines. 

We discuss the SST/IRS observations in section 2. In section 3, the data are used to 
test for a variation in the degree of ionization of the H II regions with Rq. We examine 
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the Ne/S abundance ratio for our M83 H II region sample in section 4. Section 5 describes 
how these Spitzer data are used to constrain and test the ionizing SEDs predicted by stellar 
atmosphere models. Last, we provide a summary and conclusions in section 6. 

2. Spitzer Space Telescope Observations 

In the nearly face-on (tilt 24°) spiral galaxy M83, we observed 24 H II regions, covering 
a wide range of deprojected galactocentric radii (Rg). We used the SST/IRS in the short 
wavelength, high dispersion (spectral resolution ~ 600) configuration, called the short-high 
(SH) mode (e.g., Houck et al. 2004). This covers the wavelength range from 9.9 - 19.6 /xm 
permitting cospatial observations of all four of our program emission lines: [S IV] 10.51, 
[Ne II] 12.81, [Ne III] 15.56, and [S III] 18.71 /im. The SH bandpass also covers the hydrogen 
Hua, also called H(7-6), at 12.37 /im. Unfortunately, we were not able to detect this line in 
any of the 24 H II regions observed. 

The observations were made in 2005, February and July. Figure 1 shows the regions 
and apertures observed, while Table 1 lists the H II region positions and the aperture grid 
configuration used to observe it. Nebulae with RK and deV designations are from Rumstay 
& Kaufman (1983) and de Vaucouleurs et al. (1983), respectively. The size of the SH 
aperture is 11.3"x4.7". In all cases, we chose the mapping mode with aperture grid patterns 
varying from a 1x2 grid to as large as a 2x4 grid in order to cover the bulk of the emission. 
Maps were arranged with the apertures immediately abutting each other; that is, with no 
overlap or space between them. In order to save overhead time, we clustered the objects into 
"Astronomical Observing Templates" (AOTs) with the same aperture grid pattern. 

Our data were processed and calibrated with both versions S12.0.2 and S13.2.0 of the 
standard IRS pipeline at the Spitzer Science Center (SSC). For each position, we have 16 
cycles with 30 sec ramp duration. The basic calibrated data (bed images) for each telescope 
pointing were median-combined and cleaned of rogue pixels and noisy order edges (the ends 
of the slit). The spectra were extracted and the lines measured using the IRS Spectroscopic 
Modeling, Analysis and Reduction Tool (SMART, Higdon et al. 2004). For the brighter 
lines, we find little change in the line fluxes between the two pipeline versions. Generally the 
[S IV] 10.5 line was weakest; for these we used the later version of the pipeline. The emission 
lines were measured with SMART using a Gaussian line fit. The continuum baseline was 
fit with a linear or quadratic function. Figures 2 (a)-(d) show the fits for each of the four 
lines in RK209 (object #9 in Fig. 1). A line is deemed to be detected if the flux is at least 
as large as the 3 o uncertainty. We measure the uncertainty by the product of the FWHM 
and the root-mean-square variations in the adjacent, line-free continuum; it does not include 
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systematic effects. 

We now discuss and estimate systematic uncertainties for our specific case using IRS 
in SH mode. Most likely the largest uncertainty is due to slit (aperture) loss factors. The 
pipeline flux calibration assumes that objects are point sources. Our nebulae are extended 
and that is why we mapped each with a grid that covers more than a single aperture. We did 
not make a correction for this effect. Thus we have implicitly assumed that the H II regions 
are close to the point-source limit within the SH 11.3"x4.7" aperture. If the H II region were 
uniformly extended within the SH aperture, correction factors would need to be applied to 
our fluxes. These are: 0.697, 0.663, 0.601, and 0.543 for the 10.5, 12.8, 15.6, and 18.7 /im lines 
respectively (Simpson et al. in preparation). These factors were obtained by interpolating in 
numbers provided from the blslitlossjconvert.tbl file from the Spitzer IRS Custom Extrac- 



tion tool (SPICE) for the SH module (see http:/ /ssc.spitzer.caltech.edu/postbcd/spice.html). 



For the uniformly filled aperture, the maximum uncertainty in the flux due to this effect 
would be ~46% for the [S III] 18.7 line. The correction factors would need to multiply 
our listed fluxes. We note that with regard to this effect, the fluxes listed in Table 2 are 
upper limits and that the uncertainty would be in only the direction to lower them. No 
correction factor was applied because we are likely closer to the point-source limit than the 
uniform-brightness limit. Because our science depends on line flux ratios, for our purposes, 
the possible uncertainty due to this effect would be lower, e.g., ~22% when we deal with the 
line flux ratio [S IV] 10.5/[S III] 18.7. 

According to §7.2 Spectroscopic Flux Calibration Uncertainties (Infrared Spectrograph 
Data Handbook, ver. 2.0 (http://ssc.spitzer.caltech.edu/irs/dhb), the minimum uncertainty 
possible in the absolute flux calibration of the spectroscopic products delivered by the pipeline 
is ±5% due to the "photometric uncertainty introduced by uncertainty in the angular sizes of 
the standard stars and their spectral types" . It is also stated that the maximum uncertainty 
[for a point source] is ±10% by comparisons with other Spitzer instruments. 

Any uncertainty in the flux due to a pointing error is probably small and in the worst 
case should not exceed 10%. We arrive at this estimate as follows. The absolute pointing 
accuracy of Spitzer is ~1" rms. Since a spectral map is performed by alignment of a guide 
star, this is effectively the pointing uncertainty for each spatial position in the spectral map. 
Because our smallest map dimensions are 11.3" x 9.4", an error of 1" could result in an error 
of ~10% if the source uniformly filled the map area and went to zero outside of it. 

For the brighter lines, i.e., most of the 12.8, 15.6, and 18.7 lines, the systematic 
uncertainty far exceeds the measured (statistical) uncertainty. Even for the fainter lines, we 
estimate that the systematic uncertainty exceeds the measured uncertainty. 
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In addition to the line flux, the measured FWHM and radial velocity (V r ) are listed 
in Table 2. Both the FWHM and V r are useful in judging the reliability of the line mea- 
surements. The FWHM is expected to be the instrumental width for all our lines. With 
a resolving power for the SH module of ~600, our lines should have a FWHM of roughly 
500 km s _1 . The values for V r should straddle the heliocentric systemic radial velocity for 
M83 of ~516 km s _1 (Lawrence et al. 1999). We note that the IRS pipeline does not correct 
to heliocentric radial velocities (Vheiio)- At the low ecliptic latitude of M83, the correction 
of V r to Vheiio can be almost ±30 km s _1 . It is interesting (but not pertinent) that we can 
definitely measure an average velocity differential between our February and July data of 
roughly 60 km s -1 , as expected. 



3. Variation in the degree of ionization of the H II regions with R G 

From the measured fluxes, we estimate ionic abundance ratios, including Ne ++ /Ne + , 
S 3+ /S ++ , and S ++ /Ne + , for each of the H II regions. Important advantages compared 
with prior optical studies of various other ionic ratios are: 1) the IR lines have a weak 
and similar electron temperature (T e ) dependence while the collisionally-excited optical lines 
vary exponentially with T e , and 2) the IR lines suffer far less from interstellar extinction. 
Indeed for our purposes, the differential extinction correction is negligible as the lines are 
relatively close in wavelength. In our analysis, we deal with ionic abundance ratios and 
therefore line flux ratios. In order to derive the ionic abundance ratios, we perform the 
usual semiempirical analysis assuming a constant T e and electron density (N e ) to obtain the 
volume emissivities for the four pertinent transitions. We use the atomic data described 
in Simpson et al. (2004) and Simpson et al. (in preparation) for the ions Ne + , Ne ++ , S ++ , 
and S 3+ . For the entries in Table 3, we adopt a typical value for all the M83 H II regions 
of T e = 8000 K and N e = 100 cm -3 . The H II regions in M83 are known to have a high 
metallicity (e.g., at least twice solar, Dufour et al. 1980; Bresolin & Kennicutt 2002); thus a 
value of 8000 K is not unreasonable. This is also typical of the values found for seven H II 
regions in M83 (see Tables 10 and 11 in Bresolin et al. 2005). Because of the insensitivity of 
the volume emissivities to T e , particularly when working with ratios for these IR lines, our 
results depend very little on this T e choice. The effects on our analysis due to a change in 
the assumed N e are also small as will be discussed later. 

We chose our sample of nebulae in order to cover a wide range in R G . It is a straight- 
forward geometry exercise to derive the deprojected galactocentric distances. This involves 
knowing the inclination angle (i = 24°), position angle of the line of nodes (9 = 43°), and 
distance (D = 3.7 Mpc) (de Vaucouleurs et al. 1983). We assumed the centre of the galaxy is 
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at a, 5 = 13 h 37 m 00!92, -29°51'56"7 (J2000) (source NASA/IPAC Extragalactic Database). 
Table 3 lists R G for the centre of each object. These range from 0.46 to 5.16 kpc. 

We present the variation of Ne ++ /Ne + with R G in Figure 3 using the values from 
Table 3. The error values represent the propagated flux measurement uncertainties and do 
not include the systematic uncertainties. There is extremely little change in any of these 
ratios even when using an iV e of 1000 cm" 3 , which is likely a reasonable upper limit for 
these H II regions (see Figure 9 in Bresolin et al. 2005). A linear least-squares fit indicates 
a positive correlation with R G (in kpc), 

Ne ++ /Ne+ = 0.033±0.010 + (0.011±0.0035) R G , 

with miniscule change to this equation for N e = 1000 cm -3 . For all the least-squares line fits 
in this paper, each point is given equal weight because systematic uncertainties exceed the 
flux measurement uncertainties, as discussed earlier. The positive correlation of Ne ++ /Ne + 
with R G as measured by the slope may be judged to be significant following the criterion 
that it exceeds the 3 a uncertainty. 

A similar fit to the S 3+ /S ++ vs. R G data yields 

S 3 +/S++ = -0.00079±0.0079 + (0.0067±0.0027) R G . 

There is more scatter in the plot (not shown) because the data for the [S IV] 10.5 line are 
noisier than for the other stronger lines (see Figure 2). The slope exceeds the 1 a uncertainty 
by a factor of only 2.4. Thus the increase in degree of ionization with increasing R G in this 
case would be deemed marginal. We note that one of the points (for RK268) stands out as 
far above any of the others. If we refit eliminating this far outlier, then 

S 3 +/S++ = 0.0055±0.0032 + (0.0029±0.0011) R G . 

The slope is not as steep and as 2.5 a would again be deemed marginally significant. 

Figure 4 plots the fractional ionic abundance ratio <S ++ >/<Ne + > vs. R G . This ratio 
is obtained from the S ++ /Ne + ratio by multiplying by an assumed Ne/S value (see below). 
The last three columns of Table 3 list this and other fractional ionic abundance ratios used 
in this paper. In Figure 4, the filled circles represent the points, and the solid line, the linear 
least-squares fit for an assumed N e of 100 cm" 3 . For H II regions in M83 this is likely the 
typical N e value. Here, the fit indicates a significant positive correlation with R G , 

<S ++ >/<Ne+> = 0.43±0.026 + (0.035±0.0089) R G , 

where angular brackets denote fractional ionization. The lower dashed line is the least- 
squares fit to points marked with an X derived assuming N e = 1000 cm -3 . This shows the 
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effect of higher N e on the volume emissivity of the [S III] line. The least-squares fit for this 
density is 

<S++>/<Ne+> = 0.36±0.022 + (0.030±0.0076) R G , 

and here too, the slope is statistically significant. In this figure and in these linear fits, we 
assume an Orion Nebula Ne/S abundance ratio of 14.3 (Simpson et al. 2004). Because Ne 
and S are "primary" elements, their production is expected to vary in lockstep and Ne/S 
would not be expected to show a radial gradient within a galaxy (Pagel & Edmunds 1981). 
There is a clear correlation of increasingly higher ionization with increasing R G . This is 
most likely due to the lower metallicity at larger R G causing the exciting stars to have a 
harder ionizing spectrum. The first quantitative abundance study of radial gradients for H II 
regions in M83 found dlog(0/H)/dR(kpc) = —0.09 ± 0.02 dex kpc~ x based on differential 
photoionization modeling (Dufour et al. 1980). With ~25 more years of observations, this 
gradient still appears but has flattened considerably. In their Figure 7, Bresolin & Kennicutt 
(2002) plot the O/H gradient. When we convert their units to R(kpc) using the M83 distance 
of 3.7 Mpc, we find dlog(0/H)/dR(kpc) = -0.0257 dex kpc" 1 . A shallower slope is also 
indicated in Figure 15 (see also their Figure 20) in Bresolin et al. (2005), where the M83 
H II region points are shown as open squares. In their Figure 15, M83 has the flattest O/H 
gradient among the 5 galaxies plotted. Because M83 has a bar, this tends to correlate with 
a flatter radial abundance gradient (e.g., Martin & Roy 1994) due to radial mixing. 



4. Neon to Sulfur abundance ratio 

For H II regions, we may approximate the Ne/S ratio with (Ne + + Ne ++ )/(S ++ + S 3+ ). 
This includes the dominant ionization states of these two elements. However this relation 
does not account for S + , which should be present at some level. We may safely ignore 
the negligible contributions of neutral Ne and S in the ionized region. Figure 5 shows our 
approximation for Ne/S vs. R G . There appears to be a drop in the Ne/S ratio with increasing 
R G . The expected increasing fraction of S + towards the inner galaxy regions would lead to 
a flatter gradient. Another factor that could flatten the slope is the higher dust content 
(with S, but not Ne, entering grains) expected in the inner regions due to higher metallicity 
as is the case for the Milky Way. The refractory carbonaceous and silicate grains are not 
distributed uniformly throughout the Galaxy but instead increase in density toward the 
centre. A simple model suggests the dust density is ~5 - 35 times higher in the inner parts 
of the Galaxy than in the local ISM (Sandford, Pendleton, & Allamandola 1995). The Ne/S 
abundance ratios that we derive here are considerably higher than the Orion Nebula value of 
14.3 (Simpson et al. 2004). We suspect that all our Ne/S estimates are upper limits because 



-10- 



of the two effects (not accounting for S + or dust) with the outer regions likely needing less 
of a downward correction to obtain a true Ne/S ratio. 

The premise that these estimates of Ne/S are upper limits is further supported by our 
forthcoming work with H II regions in M33. This is a similar study for which we have Spitzer 
Cycle 2 spectra of roughly 25 nebulae in the substantially face-on spiral M33. A preliminary 
report was presented at IAU Symposium 235 (Rubin et al. 2006). The Ne/S values (under 
the same approximation) fell mostly in the range from ~12 - 21. M33 H II regions have 
a lower metallicity than those in M83. Also, most of the H II regions we observed in M33 
have significantly higher ionization than those we observed in M83. These two facts tend to 
mitigate the amount of any downward correction needed to account for S + and dust. While 
there appears to be a decrease in our (approximate) Ne/S vs. R^ for M33, our data also 
indicate that the lower envelope to Ne/S is well fit by a constant value equal to the Orion 
Nebula ratio (Rubin et al. 2006 and in preparation). 



5. Constraints on the ionizing SED for the stars exciting the H II regions 

Various fractional ionic abundances are highly sensitivity to the stellar ionizing SED that 
apply to H II regions. This has been realized, for example, for the Ne + - Ne ++ ionization 
equilibrium and the total number of photons more energetic than the 41 eV Ne + ionization 
potential that are predicted by various stellar atmosphere models (e.g., Simpson et al. 2004 
and references therein). The present Spitzer data probe the Ne + and Ne ++ fractional ionic 
abundances, as well as those of S ++ and S 3+ ; thus they may be used to provide further 
constraints and tests on the ionizing SED for the stars exciting these M83 nebulae. We 
use the ratio of fractional ionizations <Ne ++ >/<S ++ > vs. <S 3+ >/<S ++ > (Figure 6a) and 
<Ne ++ >/<S 3+ > vs. <S 3+ >/<S ++ > (Figure 6b). These ionic ratios are computed using 
our photoionization code NEBULA (e.g., Simpson et al. 2004; Rodriguez & Rubin 2005). 
The lines connect the results of the nebular models calculated using the ionizing SEDs 
predicted from various stellar atmosphere models. There are no other changes to the input 
parameters, just the SED. The stellar atmospheres used are representative of several recent 
non-LTE models that apply for O-stars. We also display the results from one set of LTE 
models by Kurucz (1992). His LTE atmospheres have been extensively used in the past as 
input for H II region models; hence the comparison with the other non-LTE results reinforces 
the fact that more reliable SEDs for O-stars require a non-LTE treatment. Figures 6a and b 
dramatically illustrate how sensitive H II region model predictions of these ionic abundance 
ratios are to the ionizing SED input to nebular plasma simulations. 

For the H II region models calculated with Pauldrach et al. (2001) atmospheres, the 



- 11 - 



solid line connects models with dwarf atmospheres and the dashed line connects models with 
supergiant atmospheres. The Sternberg et al. (2003) paper also uses Pauldrach's WM-BASIC 
code. At a given T e g we have used their model with the smallest log g in order to be closest to 
the supergiant case. Because the locus using these Sternberg et al. atmosphere models is for 
the most part similar to the Pauldrach et al. supergiant locus, we do not show it in Figures 6 
to avoid clutter. For H II region models calculated with Lanz & Hubeny (2003) atmospheres 
(TLUSTY code), the solid line connects models with atmospheres with log g = 4.0 and the 
dotted and dashed lines connect models with atmospheres with log g = 3.0 to 3.5, and with 
Lyman continuum luminosities of 10 49 and 10 50 photons s -1 , respectively. The lines with 
open squares in Figures 6 are the results of our nebular models with the atmospheres in 
Martins et al. (2005) that use Hillier's CMFGEN code. 

To compare our data with the models, we need to divide the observed Ne ++ /S ++ and 
Ne ++ /S 3+ ratios by an assumed Ne/S abundance ratio. For the purposes of Figure 6, we 
adopt a constant Ne/S = 14.3, the Orion Nebula value (Simpson et al. 2004). As per 
the discussion in §4, we cannot conclude definitively whether Ne/S may vary. The open 
circles (adjusted by the assumed Ne/S) are derived from our observed line fluxes using N e 
of 100 cm~ 3 . The M83 data usually lie closest to the Pauldrach et al. supergiant loci. In 
addition the points derived from our data as well as those for the H II regions we observed in 
M33 (Rubin et al. 2006 and in preparation), which are generally of higher ionization, follow 
the trend of these theoretical loci. This is particularly notable in Figure 6b, where the other 
model loci are nearly perpendicular to the data point trend in the vicinity of where they 
intersect the data points. 

The nebular models used to generate Figures 6 are all constant density, spherical models. 
We used a constant total nucleon density (DENS) of 1000 cm" 3 that begins at the star. Each 
model used a total number of Lyman continuum photons s _1 (N Lyc ) = 10 49 . The same nebular 
elemental abundance set was used for all nebular models. We use the same "reference" set as 
in Simpson et al. (2004) because in that paper we were studying the effects of various SEDs 
on other ionic ratios and other data sets. Ten elements are included with their abundance 
by number relative to H as follows: (He, C, N, O, Ne, Si, S, Ar, Fe) with (0.100, 3.75E-4, 
1.02E-4, 6.00E-4, 1.50E-4, 2.25E-5, 1.05E-5, 3.75E-6, 4.05E-6), respectively. While 
the purpose here was not to try to match the abundances in the H II regions we observed 
in M83, the set of abundances used is roughly a factor of 1.5 higher than for Orion and not 
drastically different from solar. We have investigated the effects of changing DENS, N Lyc , 
and allowing for a central evacuated cavity, characterized by an initial radius (R^t) before 
the stellar radiation encounters nebular material. We term these shell models. In Figures 7a 
and 7b, the resulting changes to Figures 6a and 6b are shown for six nebular models run 
using two of the Pauldrach et al. (2001) supergiant atmospheres. These are listed along 
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with the 


symbol: 












(1) T cS = 


= 35000 K; 


DENS = 


1000 cm" 3 ; R init = 


= 0.5 pc; N Lyc 


i r\ A Q 

= 10 49 £ 


1 (asterisk) 


(2) T cS = 


= 35000 K; 


DENS = 


100 cm- 3 ; R init = 


pc; N Lyc = 


10 49 s- 1 


(triangle) 


(3) T eS = 


= 35000 K; 


DENS = 


100 cm" 3 ; Ri n i t = 


pc; N Lyc = 


1Q 50 s -l 


(X) 


(4) T cS = 


= 40000 K; 


DENS = 


1000 cm" 3 ; = 


= 0.5 pc; N Lyc 


= 10 49 s 


- 1 (asterisk) 


(5) T cS = 


= 40000 K; 


DENS = 


100 cm" 3 ; Rj nit = 


pc; iV L?;c = 


10 49 s- 1 


(triangle) 


(6) T cfr = 


= 40000 K; 


DENS = 


100 cm" 3 ; R^it = 


pc; N Lyc = 


10 50 s- 1 


(X) 



The original Pauldrach et al. (2001) supergiant locus and points derived from the 
Spitzer data are shown again. We also display the effects of using a different iV e to interpret 
our line measurements in terms of ionic ratios. The open stars show the results with iV e of 
1000 cm -3 while the original points (open circles) were derived with N e = 100 cm- 3 . The 
change is slight with the higher iV e shifting points to the upper right in Figure 7a and to the 
lower right in Figure 7b. 

The points for cases (3) and (6) are nearly identical to the original points for the Paul- 
drach supergiant models at the same respective T e g. This can be understood in terms of the 
ionization parameter (U), which is very useful for gauging ionization structure. An increase 
in U corresponds to higher ionization (for a given T e g). For an ionization bounded, constant 
density case, 

U = [N e N Lyc (a - a 1 ) 2 /(367rc 3 )] 1/3 , 

where (a - cki) is the recombination rate coefficient to excited levels of hydrogen, and c 
is the velocity of light (see Rubin et al. 1994, eq. 1 and adjoining discussion). Because 
(a - ai) ~ 4.10xl0- 10 T e -°- 8 cm 3 s" 1 (Rubin 1968 fit to Seaton 1959), there is only a weak 
dependence of U on T e (~ T~ - 5 ). When U is similar, as is the case here with the product 
of N e xN Lyc , the ionization structure is similar. With regard to the two shell models in 
Figures 7, the Stromgren radius is ~0.74 pc. Thus the radial thickness of the shell is slightly 
less than half the radius of the central cavity. From the visual appearance of our target H II 
regions, it is unlikely that the theoretical loci need to be tracked to higher dilutions. 

Another nebular parameter that can alter the theoretical tracks is the set of elemental 
abundances used. This is certainly well established. For instance, in the grid of model 
H II regions of Rubin (1985), the fractional ionic abundances needed here were tabulated 
considering a variation in the heavy-element abundance set of a factor of 10. The up ("U") 
and down ( "D" ) sets in that paper were meant to represent V4~0 and 1 / V4~0 times the nebular 
("N") Orion-like abundances at that time. A sampling of those models, that are closest to 
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the parameter space of interest in this paper, indicates median shifts of the following factors: 
<Ne ++ >/<S++> = 3.50, <Ne ++ >/<S 3+ > = 1.39, and <S 3+ >/<S ++ > = 2.13, when going 
from "U" to "D" sets. Much more apropos for the current situation, we have calculated two 
variants of the Pauldrach T e g- = 40000 K supergiant (canonical model in Figures 6 and 7) 
with the only change being in the nebular abundance set by scaling all the heavy elements 
by a factor of two higher and two lower than the set used for all other models. This is more 
than sufficient to cover the expected variation in our M83 sample allowing for the heavy- 
element gradient dlog(0/H)/dR(kpc) = —0.0257 dex kpc -1 (see §3). As expected, lower 
metallicity results in a shift to higher ionization. In both Figures 6a,b, the point moves 
to the upper right; the factors are <Ne ++ >/<S ++ > = 1.29, <Ne ++ >/<S 3+ > = 1.16, and 
<S 3+ >/<S ++ > = 1.11. Likewise, higher metallicity shifts the point to the lower left. This 
latter change is relatively larger with factors: <Ne ++ >/<S ++ > = 1.62, <Ne ++ >/<S 3+ > = 
1.26, and <S 3+ >/<S ++ > = 1.28. In the case of the M83 nebulae, we may conclude that 
the predicted spread in Figures 6 due to a reasonable uncertainty in nebular metallicity is 
far less than that due to the SEDs of the various stellar atmosphere models. 

There is also the effect of a change in the abundances used to compute the stellar 
atmosphere models. This will change the emergent stellar SED (e.g., Mokiem et al. 2004). As 
is the case for a change in the nebular model abundance set, such a modification in the stellar 
model will alter the shape of the SED in the same sense; that is, a higher metallicity will 
cause more opacity and soften the SED, and a lower metallicity will do the opposite. Mokiem 
et al. (2004) examined this using CMFGEN stellar models matching both the nebular and 
stellar metallicities. Their Figure 11 tracks the predicted variation in the [Ne III] 15.6/[Ne II] 
12.8 flux ratio over a range of 0.1 - 2 Z . 

Although beyond the scope of this paper, it would be interesting to compare models 
using different stellar atmosphere metallicities, especially if the environment indicates signif- 
icant departures from solar. However, presently the proper abundances are not accurately 
known and comparisons like those presented in this paper will help decipher the proper 
values. Improvements with regard to this point are deferred for a future paper, where the 
curves for different metallicities are compared with the observations. With the present paper, 
however, we have been investigating which of the models best fits the observations and such 
a comparison works most effectively with a fixed set of abundances common to all models, 
which are the solar ones. 

Such additional degrees of freedom (besides the SED) in the nebular models makes 
judging how well a particular stellar atmosphere model set fits observations more challenging. 
We will revisit these fits in a second paper that has the benefit of additional observational 
data of H II regions in M33 (Rubin et al. in preparation). 
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6. Summary and conclusions 

We have observed emission lines of [Ne II] 12.81, [Ne III] 15.56, [S III] 18.71, and [S IV] 
10.51 /iiii cospatially with the Spitzer Space Telescope using the Infrared Spectrograph (IRS) 
in short-high mode (SH). From the measured fluxes, we determined the ionic abundance 
ratios Ne ++ /Ne + , S 3+ /S ++ , and S ++ /Ne + in 24 H II regions in the substantially face-on spiral 
galaxy M83. These nebulae cover a range from 0.46 to 5.16 kpc in deprojected galactocentric 
distance Kg- We found a correlation of increasingly higher ionization with increasing Rq. 
This is seen in the variation of Ne ++ /Ne + and <S ++ >/<Ne + > with Kg (see Figures 3 and 
4). This is most likely due to the lower metallicity at larger R G causing the exciting stars to 
have a harder ionizing spectrum. 

By sampling the dominant ionization states of Ne and S for H II regions, we can ap- 
proximate the Ne/S ratio by (Ne + + Ne ++ )/(S ++ + S 3+ ). The decrease in this ratio with 
increasing R^ is more likely due to other effects than a true gradient in Ne/S. Both Ne and 
S are the products of a-chain reactions following carbon and oxygen burning in stars, with 
large production factors from core-collapse supernovae. Both are primary elements, making 
their yields depend very little on the stellar metallicity. Thus, at least to "first order" , it is 
expected that Ne/S remains relatively constant throughout a galaxy. 

As discussed in §4, our estimate for Ne/S has not accounted for the presence of S + . 
Because inner-galaxy H II regions have lower ionization (§3), it is possible that there would 
be a larger fractional ionization of S + towards the inner galaxy regions and hence a flatter 
Ne/S slope than indicated in Figure 5. A second factor which might flatten this apparent 
gradient is the higher grain content in the inner regions due to higher metallicity. All of our 
derived Ne/S abundance ratios are considerably higher than the benchmark Orion Nebula 
value of 14.3. Because we did not account for S + or dust (with the objects at larger Kg likely 
needing less of a downward correction to obtain a true Ne/S ratio), our Ne/S estimates may 
be considered upper limits. This conclusion is further supported by our subsequent work 
on H II regions in M33 where this methodology leads to lower estimates for Ne/S. This 
was briefly discussed in §4 with citations to a conference proceedings and a future paper 
(Rubin et al. 2006 and in press). Thus with observations of the set of four IR emission lines 
and the analysis we have described, there is the potential for reliable Ne/S measurements, 
especially for H II regions that have lower metallicity and higher ionization than those in 
M83. 

At the present time, the solar abundance, particularly of Ne, is the subject of much 
controversy (e.g., Drake & Testa 2005; Bahcall, Serenelli, & Basu 2006; and references in 
each of these). While we cannot directly address the solar abundance with our observa- 
tions of extragalactic H II regions, it is important to have reliable benchmarks for the neon 
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abundance. There appears to be a growing body of evidence that the Ne abundance [its 
fractional number abundance relative to hydrogen log H = 12, by definition and termed 
-A(H)] is substantially higher in the solar neighborhood, and even in the Sun itself, than 
the "canonical" solar values given in two recent, often- referenced papers. These papers have 
for the Sun: A(Ne) = 7.87, A(S) = 7.19 (Lodders 2003) and A(Ne) = 7.84, A(S) = 7.14 
(Asplund, Grevesse, & Sauval 2005). Thus Ne/S ~ 5 according to both. It is now generally 
accepted that Ne has the least well determined solar abundance among the most abundant 
elements. One of the proponents for a higher Ne abundance pointed out that an A(Ne) = 
8.29 would reconcile solar models with the helioseismological measurements (Bahcall, Basu, 
& Serenelli 2005). Using this value together with the A(S) values above, we obtain Ne/S of 
12.6 and 14.1, respectively, close to the Orion Nebula ratio (Simpson et al. 2004). 

According to calculations based on the theoretical nucleosynthesis, galactic chemical 
evolution models of Timmes, Woosley, & Weaver (1995), the Ne/S ratio in the solar neigh- 
borhood would change little, from 3.80 to 3.75, between solar birth and the present time 
(apropos for the Orion Nebula). These calculations were provided by Frank Timmes (pri- 
vate communication). With regard to these ratios being even lower than the canonical solar 
values, Timmes notes that although massive stars are expected to dominate the Ne and S 
production (and are all that were included in their non- rotating models with no wind losses), 
there is likely to be some re-distribution of Ne and S from rotation or from Wolf-Rayet phases 
of evolution, along with contributions of Ne from novae or even heavier intermediate mass 
stars. Hence, the theoretical model Ne/S ratios above should be a lower bound because some 
potential sources of Ne are missing. 

Additionally, the data here may be used as constraints on the ionizing SEDs for the stars 
exciting these nebulae by comparing the ratio of fractional ionizations <Ne ++ >/<S ++ > 
and <Ne ++ >/<S 3+ > vs. <S 3+ >/<S ++ > with predictions made from our photoionization 
models using stellar atmosphere models from several different sources. Figures 6 show the 
comparison, where we assume that the Ne/S ratio does not vary and equals the Orion Nebula 
value. Generally, the best fit is to the nebular models using the supergiant stellar atmosphere 
models (Pauldrach et al. 2001) computed with the WM-BASIC code. We note that this 
comparison is mainly qualitative since these ionic ratios depend not only on the SED, but 
also on the nebular parameters discussed as well as the effects of the stellar metallicity on the 
SED. This result is not sensitive to the electron density range, 100 and 1000 cm -3 , expected 
for these M83 nebulae. Furthermore, the points derived from our M83 data, as well as those 
for the H II regions we observed in M33 (Rubin et al. 2006 and in preparation), which are 
generally of higher ionization, follow the trend of the above-mentioned theoretical loci. In 
fact, the other model loci are nearly perpendicular to the data point trend in the vicinity of 
where they intersect the data points in Figure 6b. We reiterate that we do not infer that 
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the actual exciting stars are supergiants, but only that their SEDs have similar shape as the 
supergiant atmospheres computed using WM-BASIC. 

It is possible that a plot similar to Figures 6a,b using <Ne ++ >/<Ne + > for the ordi- 
nate will also have utility in comparing data with models using different SEDs. While the 
<Ne ++ >/<Ne + > ratio has the advantage of being independent of elemental abundance ra- 
tios, it appears to be more sensitive to the nebular parameters than does the <Ne ++ >/ <S 3+ > 
ratio. This fact tends to make it less unique in its ability to discriminate between the stellar 
SEDs we present in this paper. We defer this investigation until our next paper that has the 
addition of more Spitzer data with H II regions in M33. 

This work is based on observations made with the Spitzer Space Telescope, which is 
operated by the Jet Propulsion Laboratory, California Institute of Technology under NASA 
contract 1407. Support for this work was provided by NASA for this Spitzer program 
identification 3412. RHR had further support from the NASA Long-Term Space Astrophysics 
(LTSA) program. We thank Frank Timmes for providing information on the Ne/S ratio from 
a nucleosynthesis, galactic chemical evolution perspective. We thank Brendan Wakefield and 
Danny Key for assistance with the data reduction. Valuable comments by the referee are 
much appreciated. The nebular models were run on a Cray computer at JPL. Funding for 
its use in this investigation was provided by the JPL Office of the Chief Information Officer. 
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Fig. 1. — The positions and apertures observed for 24 H II regions are shown in red super- 
imposed on an Ha image of the nearly face-on (tilt 24°) M83. The nebulae are numbered W 
to E (see Table 1). N is up and E left. 
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Fig. 2. — Measurements of the four emission lines in the H II region RK209 (#9 in Fig. 1): 
(a) [S IV] 10.5 /im; (b) [Ne II] 12.8 /im; (c) [Ne III] 15.6 /im; and (d) [S III] 18.7 fim. The 
data points are the filled circles. The fits to the continuum and Gaussian profiles are the 
solid lines. Such measurements provide the set of line fluxes for further analysis. 
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1 2 3 4 5 6 

R G (deprojected) (kpc) 

Fig. 3. — Plot of the ionic abundance ratio Ne ++ /Ne + , which is derived from the measured 
line flux ratios for each H II region, vs. R^. We assume an electron density (N e ) of 100 cm -3 . 
There is extremely little change with N e over the range expected for these regions (see text). 
The linear least-squares fit indicates a positive correlation with R^. Error bars here and in 
Figures 4 and 5 are for the propagated measurement uncertainties and do not include the 
systematic uncertainties (see text). 



-23- 




Fig. 4. — Plot showing the fractional ionic abundance ratio <S ++ >/<Ne + > vs. R G . Circles 
represent the points and the solid line the linear least-squares fit for an assumed N e of 
100 cm -3 . For H II regions in M83 this is likely the typical N e value with an upper limit 
of roughly 1000 cm~ 3 . We plot with an X the points and a dashed line the least-squares 
fit for an assumed N e = 1000 cm" 3 , thereby showing the effect of higher N e on the volume 
emissivity of the [S III] line. The plotted <S ++ >/<Ne + > ratio assumes an Orion Nebula 
Ne/S abundance ratio of 14.3 (Simpson et al. 2004). Because Ne and S are "primary" 
elements, their production is expected to vary in lockstep and Ne/S would not be expected 
to show a radial gradient within a galaxy (Pagel & Edmunds 1981). There is less scatter 
about the fit than for the Ne ++ /Ne + ratio because the [Ne III] line is weaker than either 
the [Ne II] or [S III] line in these objects. Error bars are shown for the solid points only to 
avoid congestion. 





1 2 3 4 5 6 

R G (deprojected) (kpc) 



Fig. 5.— Ne/S, as approximated by (Ne+ + Ne ++ )/(S ++ + S 3+ ) (see text) vs. R G . These 
ratios are all larger than the Orion Nebula value of 14.3. There also appears to be a drop in 
the Ne/S ratio with increasing R^. However we do not account for S + , which must exist in 
these low ionization M83 H II region. Accounting for S + would lower our estimate of Ne/S. 
The expected increasing fraction of S + towards the inner galaxy regions would also lead to 
a flatter gradient. Another factor that could lower all the derived Ne/S values, as well as 
flatten the slope, is the higher dust content (with S, but not Ne, entering grains) expected 
in the inner regions due to higher metallicity. 
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Fig. 6. — a) Theoretical predictions of the fractional ionization ratios <Ne ++ >/<S ++ > vs. 
<S 3+ >/<S ++ >, computed using our photoionization code NEBULA. The lines connect the 
results of nebular models calculated with the ionizing SEDs predicted from various stellar 
atmosphere models as labeled, changing no other parameter except the SEP . For the H II 
region models calculated with Pauldrach et al. (2001) atmospheres, the solid line connects 
models with dwarf atmospheres and the dashed line connects models with supergiant at- 
mospheres. For the H II region models calculated with Lanz & Hubeny atmospheres, the 
solid line connects models with atmospheres with log g = 4.0 and the dotted and dashed 
lines connect models with atmospheres with log g = 3.0 to 3.5, and with Lyman contin- 
uum luminosities of 10 49 and 10 50 photons s _1 , respectively. To compare our data with 
the models, we need to divide the observed Ne ++ /S ++ and Ne ++ /S 3+ ratios by an assumed 
Ne/S abundance ratio. We use the Orion Nebula Ne/S = 14.3. The open circles (adjusted 
by the assumed Ne/S) are derived from our observed line fluxes using iV e of 100 cm -3 , 
b) The same as panel a) except the ordinate is <Ne ++ >/<S 3+ >. Both panels dramatically 
illustrate the sensitivity of the H II region model predictions of these ionic abundance ratios 
to the ionizing SED that is input to nebular plasma simulations. The M83 data, for the 
most part, appear to lie closest to the Pauldrach et al. supergiant loci. 
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Fig. 7. — a) This is similar to Figure 6a. We again show the locus that uses the Paul- 
drach et al. (2001) supergiant atmospheres. Here we display the results of making some 
changes to the nebular parameters for the T cS = 35000 and 40000 K stellar atmospheres. 
The points with an * are for a model with a central cavity of radius 0.5 pc (see text); 
those with a triangle have a density of 100 instead of 1000 cm~ 3 ; and those with an X 
have a larger number of Lyman continuum photons s _1 (Nl vc ) = 10 50 instead of 10 49 . 
Points corresponding to the same T c g are connected by the solid lines. Again, the open 
circles are derived from our observed line fluxes using N e of 100 cm -3 . The open stars 
show the results of using N e = 1000 cm -3 - a slight shift of the points to the upper right, 
b) This is similar to Figure 6b with the same modifications as described for panel a). As 
shown by the open stars, an N e = 1000 cm" 3 (compared with N e = 100 cm -3 ) shifts the 
points slightly to the lower right. 
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Table I. H II Regions Observed in M83 



Order 


H II Region 


RA J2000 


DEC 


Aperture Grid 


1 


RK275 


13 36 40.3 


-29 51 21 


1x2 


2 


RK268 


13 36 42.5 


-29 52 09 


1x2 


3 


RK266 


13 36 43.4 


-29 52 22 


2x2 


4 


RK230 


13 36 50.7 


-29 52 02 


1x3 


5 


deVIO 


13 36 52.3 


-29 53 00 


2x3 


6 


RK213 


13 36 52.6 


-29 51 46 


1x3 


7 


deV13 


13 36 52.7 


-29 52 46 


2x3 


8 


RK211 


13 36 52.9 


-29 51 11 


2x4 


9 


RK209 


13 36 53.2 


-29 51 31 


2x4 


10 


RK201 


13 36 53.9 


-29 48 54 


2x4 


11 


RK198 


13 36 54.3 


-29 50 47 


2x4 


12 


deV22 


13 36 54.7 


-29 53 05 


2x4 


13 


RK172 


13 36 57.1 


-29 51 55 


2x2 


14 


RK154 


13 36 58.7 


-29 48 06 


1x2 


15 


deV28 


13 36 59.0 


-29 51 26 


1x3 


16 


deV31 


13 37 00.0 


-29 52 19 


1x3 


17 


RK137 


13 37 01.4 


-29 51 27 


2x4 


18 


RK135 


13 37 02.0 


-29 55 31 


2x2 


19 


RK120 


13 37 03.5 


-29 54 02 


1x3 


20 


RK110 


13 37 04.7 


-29 50 58 


2x3 


21 


RK86 


13 37 07.1 


-29 49 36 


1x3 


22 


RK69 


13 37 08.5 


-29 52 04 


1x3 


23 


deV52+RK70 


13 37 08.6 


-29 52 11 


1x4 


24 


RK20 


13 37 16.9 


-29 53 14 


2x3 
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Table 2. M83 Line Measurements 



Order 


Source 


Line 


Flux 


la error 


FWHM 


V r 






^iin 


W cm- 2 


W cm- 2 


km s -1 


km s -1 


1 


RK275 


10.5 


9.73E-22 


1.19E-22 


360 


454 






12.8 


2.25E-20 


1.96E-22 


510 


627 






15.6 


3.27E-21 


1.12E-22 


493 


680 






18.7 


1 A A T7^ o r\ 

1.44E-20 


1.56E-22 


501 


610 


2 


RK268 


10.5 


1.15E-21 


2.16E-22 


809 


613 






12.8 


4.06E-21 


1.60E-22 


504 


651 






15.6 


1.14E-21 


4.15E-23 


456 


693 






18.7 


2.41E-21 


9.71E-23 


465 


606 


3 


RK266 


10.5 


1.62E-21 


4.20E-22 


711 


650 






12.8 


2.64E-20 


3.12E-22 


504 


657 






15.6 


5.22E-21 


5.90E-22 


502 


690 






18.7 


1.73E-20 


3.65E-22 


496 


631 


4 


RK230 


10.5 


8.14E-22 


1.22E-22 


385 


417 






12.8 


1.65E-20 


2.47E-22 


526 


638 






15.6 


4.23E-21 


1.84E-22 


513 


694 






18.7 


9.61E-21 


3.23E-21 


505 


629 


5 


deVIO 


10.5 


2.13E-21 


1.22E-22 


290 


467 






12.8 


8.41E-20 


6.42E-22 


483 


674 






15.6 


1.02E-20 


1.15E-21 


517 


718 






18.7 


4.30E-20 


1 A O TV O 1 

1.43E-21 


498 


653 


6 


RK213 


10.5 


7.40E-22 


1.23E-22 


392 


508 






12.8 


8.51E-20 


2.78E-22 


476 


654 






15.6 


5.68E-21 


2.72E-22 


479 


697 






18.7 


4.59E-20 


3.73E-22 


487 


635 


7 


deV13 


10.5 


2.01E-21 


2.66E-22 


293 


461 






12.8 


9.80E-20 


7.14E-22 


472 


672 






15.6 


1.05E-20 


8.54E-22 


502 


706 






18.7 


5.21E-20 


4.77E-22 


491 


648 


8 


RK211 


10.5 


2.30E-21 


2.25E-22 


326 


434 






12.8 


1.01E-19 


1.15E-21 


499 


624 






15.6 


1.03E-20 


4.90E-22 


529 


665 






18.7 


5.95E-20 


2.32E-21 


495 


605 


9 


RK209 


10.5 


2.90E-21 


3.25E-22 


404 


403 






12.8 


1.68E-19 


1.04E-21 


467 


580 






15.6 


1.95E-20 


4.78E-22 


491 


631 






18.7 


9.77E-20 


5.97E-22 


478 


563 


10 


RK201 


10.5 


2.79E-21 


4.19E-22 


495 


434 






12.8 


2.76E-20 


4.30E-22 


484 


566 






15.6 


5.73E-21 


5.53E-22 


906 


587 






18.7 


1.58E-20 


1.53E-21 


486 


567 


11 


RK198 


10.5 


3.45E-21 


3.99E-22 


548 


362 






12.8 


7.25E-20 


1.02E-21 


479 


549 






15.6 


1.22E-20 


3.01E-22 


516 


587 






18.7 


4.71E-20 


8.02E-22 


487 


528 


12 


deV22 


10.5 


2.99E-21 


4.76E-22 


716 


280 
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Table 2 — Continued 



Order 


Source 


Line 
fim 


Flux 

W cm^ 2 


la error 
W cm' 2 


FWHM 
km s — 1 


Mr 

km s -1 






12.8 


2.00E-19 


1.91E-21 


490 


512 






15.6 


1.77E-20 


3.91E-22 


469 


655 






18.7 


9.60E-20 


9.62E-22 


491 


595 


13 


RK172 


10.5 


4.67E-22 


1.20E-22 


343 


539 






12.8 


1.03E-20 


3.00E-22 


514 


605 






15.6 


1.57E-21 


2.57E-22 


431 


658 






18.7 


4.59E-21 


2.21E-22 


437 


601 


14 


RK154 


10.5 


6.78E-22 


9.92E-23 


300 


385 






12.8 


1.11E-20 


1.50E-22 


502 


551 






15.6 


2.52E-21 


1.30E-22 


550 


570 






18.7 


7.29E-21 


2.2bhj-22 


506 


530 


15 


dcV28 


10.5 


5.21E-22 


1.51E-22 


497 


407 






12.8 


1.39E-20 


2.81E-22 


488 


563 






15.6 


2.01E-21 


2.31E-22 


809 


593 






18.7 


6.44E-21 


2.82E-22 


470 


560 


16 


deV31 


10.5 


1.04E-21 


7.98E-23 


294 


463 






12.8 


1.14E-19 


6.59E-22 


496 


621 






15.6 


8.49E-21 


4.04E-22 


542 


661 






18.7 


4.19E-20 


O O A T7"> no 

2.84E-22 


505 


602 


17 


RK137 


10.5 


1.56E-21 


O / • ^ Tr* no 

3.67E-22 


415 


342 






12.8 


1.71E-19 


1.09E-21 


488 


518 






15.6 


1.30E-20 


5.88E-22 


557 


570 






18.7 


9.05E-20 


7.71E-22 


493 


502 


18 


RK135 


10.5 


7.76E-22 


1.15E-22 


325 


531 






12.8 


2.38E-20 


1.04E-21 


495 


649 






15.6 


1.95E-21 


2.40E-22 


466 


661 






18.7 


1.42E-20 


1.01E-21 


480 


629 


19 


RK120 


10.5 


1.99E-21 


8.66E-23 


439 


419 






12.8 


5.66E-20 


/I A A T7"> no 

4.44E-22 


484 


604 






15.6 


6.84E-21 


1.72E-22 


474 


646 






18.7 


3.23E-20 


3.00E-22 


493 


599 


20 


RK110 


10.5 


1.43E-21 


3.02E-22 


458 


379 






12.8 


1.13E-19 


9.78E-22 


480 


499 






15.6 


9.85E-21 


3.98E-22 


523 


526 






18.7 


4.86E-20 


4.56E-22 


476 


483 


21 


RK86 


10.5 


2.06E-21 


1.51E-22 


500 


337 






12.8 


1.03E-19 


6.76E-22 


498 


486 






15.6 


9.74E-21 


1.44E-22 


486 


483 






18.7 


5.66E-20 


3.38E-22 


493 


451 


22 


RK69 


10.5 


8.41E-22 


1.89E-22 


582 


226 






12.8 


3.32E-20 


2.97E-22 


467 


534 






15.6 


4.18E-21 


2.22E-22 


516 


545 






18.7 


1.54E-20 


3.84E-22 


477 


498 


23 


deV52+RK70 


10.5 


9.80E-22 


1.44E-22 


352 


421 






12.8 


7.08E-20 


8.77E-22 


487 


518 
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Table 2 — Continued 



Order Source 


Line 


Flux 


lcr error 


FWHM 


V r 




^iin 


W cm- 2 


W cm" 2 


km s -1 


km s -1 




15.6 


7.03E-21 


2.78E-22 


541 


555 




18.7 


3.61E-20 


3.64E-22 


496 


489 


24 RK20 


10.5 


6.61E-22 


1.60E-22 


328 


371 




12.8 


1.78E-20 


4.29E-22 


485 


512 




15.6 


3.87E-21 


3.19E-22 


580 


538 




18.7 


7.85E-21 


4.53E-22 


450 


499 



Table 3. Derived Parameters for the H II Regions in M83 



Order 


Source 


kpc 


Ne+ 
S++ 




Ne++ 
S++ 




Ne++ 
Ne+ 




g3+ 

S++ 






Ne 
S 






<s++> 

<Ne+> 


<Nc++> 
<S++> 


<Nc++> 
<S 3+ > 


i 
l 


rxivz f 


k 1 ft 
0. id 


OQ A 1 n 
zo. 4±U. 


Q 
O 


1.4y±U 


UO 


n nftQ K_i_n 
U.UooOitU 


nnnq 
.UUzo 


n ni a n 1 n 
U.U14o±U 


nni 

.UUlo 


-1 
z4 


ft 1 n 


Q 
,0 


U 


fti n_Ln nnQ 
.DlUztU.UUo 


n 1 c\a J_n nn/i 
U. 1U4±U.UU4 


7 oftj_n no 

/ .zDztu.yz 


o 

Z 


T3 Ti'OftQ 
mVZDo 


A -1 
4.01 


OK A -L 1 
ZO.4r.tl. 


A 
4 


q 1 n~un 
0. lUztU 


1 7 
1 1 


n 1 oo-Un 
U. IzzztU 


.UU ( 


n im ~un 

U. lUlltU 


ni 
uiy 


ZD 


. o±l 


U 


n 
U 


Kft/1-Un nQi 
.0D4±U.Uol 


n 01 7-un ni 
U.zl / ztU.Ulz 


O 1 K-Un A 1 
Z. 10ltU.41 


Q 
O 


T3 Iif Oftft 


A OQ 

4.zy 


zz.yztu. 


K 



1 no_i_n 

i.yoztu. 


OQ 

zo 


U.UoOOltU 


.uuyo 


n ni oo-Un 
u.uiyy±u 


nnKO 
UUOz 


OA 
Z4 


a ~un 


,0 


U 


ftOQ-un ni k 
.OzoitU.UlO 


n 1 QO-un ni ft 
U. loy±U.UlD 


ft Qft_|_1 Oft 

D.yo±i.yo 


A 

4 


ItlYzoU 


Kn 

z.ou 


ok <7 1 n 
zO. / ±U. 


n 

y 


001 n 
z.ooztU. 


1 ft 
10 


n 1 1 i n 
U. llzztU 


nne 

.UUO 


n ni on 1 n 
U.UloU±U 


nnoo 
UUzo 


ZO. 


1 i n 


n 

,y 


U 


KKftj_n non 
.OODitU. UzU 


n ono_Ln ni 1 
U.zUzzLU.Ull 


1 1 OJ-1 Q 

ll.zztl.o 


K 



cie v iu 


Q/l 
Z.o4 


OQ 7-U1 
ZO. f ±1. 


A 

u 


1 KK-Un 

l.OoztU. 


1 Q 
. lo 


n nK/in-un 
U.U04U±U 


nnfti 

.UUDl 


n ni nK-un 
U.UlUOztU 


nnnv 
UUU 1 


zy 


Q_|_1 


n 

U 


n 

U 


a OQ-Un ni ft 
.4yo±U.UlD 


n 1 nQ-un ni q 
U. lUoztU.Ulo 


1 n Q-U1 Q 

1U. ozt l.o 


ft 
O 


"D I^OI Q 
ItlVzlo 


z.U / 


07 n 1 n 
z ( .y±u. 


n 
z 


n qi q 1 n 
U.oloitU 


uoy 


n nonoxn 
U.UzyzztU 


nni a 
.UU14 


n nnQ/i Qj_n 
U.UUo4o±U 


nnnKV 
UUUO f 


OQ 
ZO. 


ft 1 n 



z 


n 
U 


ki Qj_n nn/i 
.0lo±U.UU4 


n nKftnj_n nnoc 
U.UODy±U.UUzo 


1 ft ftj_o n 

lo.oztz.y 


>7 
I 


ele v lo 


OIK 

z. 10 


oq 1 n 
zo.zztU. 


Q 






.11 


U.U4oy±U 


nnoo 
UUoo 


n nnoi i n 
U.UUolo±U 


nni nn 

uuiuy 


on 

zy 


1 n 
.o±U 


Q 
O 


U 


Knftj_n nnft 
.OUDztU.UUo 


n nnoK_i_n nn7ft 
U.UyzOztU.UU ( D 


1 1 Q-i-1 Q 

ll.o±l.o 


Q 

O 


IxlYzl 1 


Z.Zo 


OK /I -U1 
Z0.4±l. 


u 


l.loztU. 


C\1 
\) ( 


U.U440±U 


nnoo 

.UUzZ 


n nncon-un 
U.UUozUztU 


nnncfi 
UUUoD 


Oft 

zo 


Q_|_1 

.ozt 1 


n 

U 


n 
U 


KftQ-un noo 
.ODoztU.Uzz 


n n7004-fi nn/iQ 
u.u / yz±u.uu4o 


O fift_l_1 OK 

y.oDiti.uo 


n 

y 


ItlYzUy 


no 


ok 1 n 
zO.ozbU. 




z 


1 Q1 I n 


no 

Uo 


n nKn7_i_n 
U.UOU f ±U 


nni q 

UUlo 


n nnftQn_i_n 


nnn7i 
UUU ( 1 


07 

Z ( 


n 1 n 
.U±U 



,z 


n 
U 


KKQj_n nnK 
.OOo±U.UUO 


n nm 7J_n nnoQ 
U.Uyi f ±U.UUzv5 


1 A ft_L1 T 
14.0±1. ( 


1U 


T5 Lf om 
xvlYzUl 


4.UU 


Oft /I -UO 
Z0.4ztZ. 


K 



O QO-Ufl 

z.oy±u. 


Q 


u.uyu / ±Vj 


nnco 
.UUoy 


n nQ7ft~un 
U.Uo 1 DztU 


nnftv 
UUD 1 


7 
Z ( . 




ft 
.0 


U 


K/i O-Un nKO 
.04z±U.UOz 


n 1 ft7_un noQ 
U. ID / ztU.Uzo 


A A K-UO 70 

4.4o±u. i y 


1 "1 


T5 LT1 OQ 


z. lo 


oq 1 -Un 
zo. litu. 


K 



l. / UztU 


nK 
UO 


n nvQQ-Ufi 
U.U / ooztU 


nnoi 
.UUzl 


n ni KK-Un 
U.UlOOztU 


nni s 

UUlo 


OA 
Z4 


.OztU 


.4 


n 

U 


fti Q-Un ni q 

.DloztU.Ulo 


n 1 1 Q-Un nn/i 
U.liy±U.UU4 


7 fi/i_i_n on 
1 .D4±u.yu 


1 

Iz 


TOO 

ue v zz 


1 01 

i.yi 


q 1 a i n 
ol.4±U. 


1 


1 01 i n 
l.zlztU 


no 
Uo 


U.Uoo f ±U 


nnnn 

uuuy 


n nnftfto_Ln 


nni nK 
UU1UO 


QO 
OZ 


a 1 n 
4±U 


Q 
O 


n 
U 


a Kft_i_n nnft 
.40D±U.UUD 


n OQ/iQ-Ln nnon 
U.Uo4y±U.UUzU 


1 O Q-LO 1 

Iz.oiLz.l 


13 


RK172 


0.95 


32.1±1. 


.8 


2.31±0. 


.38 


0.0719±0 


.0117 


0.0216±0 


.0057 


33 


.7±1. 


,6 





.445±0.025 


0.162±0.027 


7.49±2.27 


14 


RK154 


4.41 


22.9±0. 


.8 


2.27±0. 


.14 


0.0992±0 


.0053 


0.0198±0 


.0030 


24 


,7±0. 


,8 





.623±0.021 


0.159±0.009 


8.05±1.25 


15 


deV28 


0.77 


32.5±1. 


.5 


2.05±0 


25 


0.0630±0 


.0074 


0.0172±0 


.0050 


34 


.Oil 


.5 





.439±0.021 


0.143±0.018 


8.35±2.60 


16 


dcV31 


0.46 


40.8±0. 


.4 


1.33±0. 


06 


0.0326±0 


.0016 


0.00525±0. 


.00041 


41 


,9±0. 


,3 





.350±0.003 


0.0931±0.0045 


17.7±1.6 


17 


RK137 


0.57 


28.4±0. 


3 


0.943±0. 


043 


0.0332±0 


.0015 


0.00366±0. 


00086 


29 


3±0. 


,2 





.503±0.005 


0.0660±0.0030 


18.0±4.3 


18 


RK135 


4.05 


25.1±2. 


1 


0.898±0 


.127 


0.0357±0 


.0047 


0.0116±0 


.0019 


25 


.7±1. 


.8 





.568±0.046 


0.0629±0.0089 


5.43±1.05 


19 


RK120 


2.48 


26.3±0. 


3 


1.39±0. 


.04 


0.0528±0 


.0014 


0.0131±0 


.0006 


27 


.3±0. 


3 





.543±0.006 


0.0973±0.0026 


7.43±0.37 


20 


RK110 


1.38 


34.9±0. 


.4 


1.33±0. 


.06 


0.0381±0 


.0016 


0.00624±0. 


.00132 


36 


.0±0. 


,3 





.409±0.005 


0.0931±0.0039 


14.9±3.2 


21 


RK86 


2.93 


27.3±0. 


2 


1.13±0 


02 


0.0413±0 


.0007 


0.00774=b0 


.00057 


28 


.3±0. 


.2 





.523±0.005 


0.0791±0.0013 


10.2±0.8 


22 


RK69 


1.89 


32.5±0. 


8 


1.78±0. 


.10 


0.0550±0 


.0030 


0.0116±0 


.0026 


33 


.8±0 


.8 





.440±0.011 


0.125±0.007 


10.8±2.5 


23 


dcV52+RK70 


1.92 


29.4±0. 


.5 


1.28±0. 


05 


0.0434±0 


.0018 


0.00576±0. 


00085 


30 


.5±0. 


3 





.486±0.008 


0.0894±0.0036 


15.5±2.4 


24 


RK20 


4.28 


34.0±2. 


1 


3.23±0 


32 


0.0950±0 


.0082 


0.0179±0 


.0044 


36 


.6±2 


.1 





.420±0.026 


0.226±0.023 


12.7±3.2 



